Asymmetric dark matter and the Sun 
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Cold dark matter particles with an intrinsic matter-antimatter asymmetry do not annihilate 
after gravitational capture by the Sun and can affect its interior structure. The rate of capture 
is exponentially enhanced when such particles have self-interactions of the right order to explain 
structure formation on galactic scales. A 'dark baryon' of mass 5 GeV is a natural candidate and 
has the required relic abundance if its asymmetry is similar to that of ordinary baryons. We show 
that such particles can solve the 'solar composition problem'. The predicted small decrease in the 
low energy neutrino fluxes may be measurable by the Borexino and SNO+ experiments. 
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We consider the capture by the Sun of asymmetric dark 
matter (ADM) particles which have a relic asymmetry 
just as do baryons, in contrast to the usual candidates 
for cold dark matter (CDM) such as supersymmetric 
neutralinos which have a relic thermal abundance deter- 
mined by 'freeze-out' from chemical equilibrium. Hence 
ADM does not annihilate upon capture in astrophysical 
bodies such as the Sun, leading to a build up of its con- 
centration. In particular, self-interactions can lead to an 
exponential increase of the ADM abundance in the Sun 
as it orbits around the Galaxy, accreting dark matter. 

ADM does not have the usual indirect signatures e.g. 
there will be no high energy neutrino signal from annihi- 
lations in the Sun. Instead ADM will alter heat transport 
in the solar interior thus affecting the low energy neutrino 
flux. This had been proposed as a solution to the 'solar 
neutrino problem' [1-.--3J. Although the solution is now 
understood to be neutrino oscillations [I], small changes 
induced by accreted CDM particles may account for the 
current discrepancy [5] between helioseismological data 
and the revised 'Standard Solar Model' (SSM). 

An asymmetry in ADM similar to that in baryons nat- 
urally explains why their observed abundances are of the 
same order of magnitude. If ADM arises from a strongly 
coupled theory (like the baryon of QCD), then there is a 
conserved U(l) global symmetry (like B number in QCD) 
which guarantees stability of the lightest U(l) charged 
object. Technicolour models of electroweak symmetry 
breaking [6\ provide an example of ADM from strong dy- 
namics in the form of the lightest neutral technibaryon 
[7] . Recently, new viable types of technibaryon dark mat- 
ter (TIMPs) [S] as well as other particle candidates for 
ADM [5] have been suggested. 

While gravitational instability in collisionless CDM 
provides a good explanation for the large-scale struc- 
tures of the universe, observations on galactic and smaller 
scales suggest that CDM may be self-interacting jTI] [T2] . 
If ADM arises as a composite TIMP or a 'dark baryon' 
from a strongly coupled model, then it would naturally 
have such self-interactions. 

We consider how the capture of self-interacting ADM 
by the Sun can alter helioseismology and low energy neu- 
trino fluxes. 



CAPTURE OF SELF-INTERACTING ADM 

We refer to earlier discussions of the capture by the 
Sun of heavy Dirac neutrinos having an asymmetry [13] , 
and of symmetric CDM with self-interactions [T4 . The 
capture rate for CDM particles \ with both an asymmetry 
and self-interactions is governed by the equation: 
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Here C x j^ is the usual rate of capture of CDM particles 
by scattering off nuclei (mainly protons) within the Sun, 
while C xx is the rate of self-capture through scattering 
off already captured \ particles. Hence the number of 
captured particles would have grown as 
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cross-section for self-captures cannot increase beyond nr 



i.e. exponentially for t > C * . However the effective 
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where r x is the scale-height of the region where they are 
gravitationally trapped pp. The linear growth by con- 
trast can continue up to the 'black disk' limit i.e. ttRq. 
In both cases there is an additional enhancement due to 
'gravitational focussing' [TJ [T3] as we quantify later. The 
ejection of captured ADM particles by recoil effects in 
the self-scattering can be neglected [14] and evaporation 
is negligible for a mass exceeding 3.7 GeV [15] , 

The ADM capture rate through spin-independent (SI) 
and spin-dependent (SD) interactions can be written [16) : 
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where cf M = 6.4x lO^s" 1 and = i.7 x 10 25 s -1 , piocai 
is the estimated local CDM density, and FiimS) encodes 
the form factors for different nuclei i weighted by the 
solar chemical composition — the sum is over all nuclei 
(hydrogen only for SD interactions). Here c^Af is the 
ADM- nucleus cross-section, which is related to ct x n, the 
ADM-nucleon cross-section [16 . For spin-independent 
interactions, ct x n is constrained by direct detection ex- 
periments such as CDMS-II Q7], XENON10 [TS] and Co- 
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GeNT jUj to be < 10~ 39 cm 2 for m x = 5 GeV. For spin- 
dependent interactions the constraints are considerably 
weaker, e.g. PICASSO [2U] sets the strongest bound of 

< IQ-36 cm 2 f or i-jjjg mass> 

Next we consider the self-capture rate in the Sun [13] : 

s ■ ^csc 

(R e ) ... erf (tj) 
xx = \l2 Plocal x 5 ^ ^ 

where s x = a xx /m x is the ADM self-interaction cross 
section divided by its mass, and t> csc (i?o) ~ 618 krns" 1 
is the escape velocity at the surface of the Sun, which 
is assumed to be moving at w = 220 kms -1 through 
a Maxwcll-Bolzmann distribution of CDM particles with 
velocity dispersion, v ~ 270 kms -1 . Here ((f)) ~ 5.1 is the 
average over tf>(r) = v 2 sc (r) j 'v 2 sc (R) and r\ = x/3/2vq/v. 

Self-interacting CDM was proposed [TT] to account 
for observations of galactic and subgalactic structure on 
scales < a few Mpc which are not in accord with numer- 
ical simulations using collisionless cold particles. The 
discrepancy can be solved if CDM has a mean free path 
against self- interactions of A ~ 1 kpc— 1 Mpc correspond- 
ing to a self-scattering cross-section between s x ~ 8 x 
10" 22 and - 8 x 10~ 25 cm 2 GeV _1 [TT] . A detailed anal- 
ysis sets an upper limit of s x < 10~ 23 cm 2 GeV^" 1 [12], 
while a study [21] of the colliding 'Bullet cluster' of galax- 
ies implies a stronger bound of ~ 2 x 10~ 24 cm 2 GeV _1 , 
which we adopt for our calculations below. 

A 'dark baryon' from a QCD-like strongly interacting 
sector but with a mass of about 5 GeV is a natural candi- 
date for ADM. Its relic density is linked to the relic den- 
sity of baryons via fl x ~ (rn x M x /m^NB)^B where Nb,x 
are the respective asymmetries. If Mb ~ M x (e.g. if both 
asymmetries are created by 'leptogenesis' [22]) then the 
required CDM abundance is realised naturally. The self- 
interaction cross-section of such a neutral particle can 
be estimated by scaling up the neutron self-scattering 
cross-section ~ 10~ 23 cm 2 [23] as: <j xx = {m n /m x ) 2 a nn 
which is just of the required order. Note that the self- 
annihilation cross-section will be of the same order which 
ensures that the ADM thermal (symmetric) relic abun- 
dance is negligible, just as it is for baryons. 

Photon exchange, via a magnetic moment of the dark 
baryon, will give rise to both spin-independent and spin- 
dependent interactions of x with nucleons. Recently this 
has been investigated in a model of a 5 GeV dark baryon 
in a 'hidden sector' interacting with the photon through 
mixing with a hidden photon magnetic moment [10] . 
From this model we infer that spin-independent cross- 
section with nuclei of 0(1O~ 39 ) cm 2 can be achieved. 
Moreover this will be accompanied by spin-dependent in- 
teractions which would aid further in the heat transport 
in the Sun as discussed below. Since the photon couples 
only to the proton in direct detection experiments, the 
limit on <t x n is degraded for this model to ~ 4x 10~ 39 cm 2 
which we adopt as an example later. 



HELIOSEISMOLOGY AND SOLAR NEUTRINOS 

Fig. [T] shows the growth of the number of captured 
ADM particles in ratio to the number of baryons in the 
Sun, for a scattering cross-section on nucleons as large 
as is experimentally allowed, including the 'gravitational 
focussing' factor of (v csc (r)/v) 2 [15] and setting r = R Q 
or r x (~ O.O7i? for m x = 5 GeV) as appropriate. 




Age of Sun [Gyr] 

FIG. 1: Growth of the relative abundance of 5 GeV mass 
ADM particles in the Sun until its present age (vertical line) 
assuming s xx = 2 x 10 _24 cm 2 GeV _1 , and a xN = 10 -39 cm 2 
(solid line) and 10~ 36 cm 2 (dashed line), these being the maxi- 
mum experimentally allowed values for spin-independent and 
spin-dependent interactions respectively. Also shown is the 
'black disk' limit (dotted line) for the Sun. 

Note that due to the self-captures, the limiting abun- 
dance N x /Nq ~ 2 x 10~ n is almost independent of the 
actual scattering cross-section. Such an ADM fraction in 
the Sun can affect the thermal conductivity and thereby 
solar neutrino fluxes [lj [2] . The SSM 24J predicts 3 times 
the observed neutrino flux (the 'Solar neutrino problem') 
but this is now well understood taking into account neu- 
trino oscillations [3J. Moreover until recently, the SSM 
with the 'standard' solar composition [25] agreed very 
well with helioseismology [26] . However the revision of 
the solar composition [27] means that the SSM no longer 
reproduces the sound speed and density profile so there 
is now a 'solar composition problem' [5]. We show that 
the presence of ADM in the Sun can resolve this problem 
and precision measurements of solar neutrino fluxes can 
constrain the properties of self-interacting ADM. 

A simple scaling argument gives for the luminosity car- 
ried by the ADM [TJ: 

£ x ~4xlO 12 L ^^i /^, (5 ) 
x °N Q a Q V m x ' 

where Lq ~ 4 x 10 33 ergss _1 . When the ADM mean free 
path X x is large compared to the scale-height r x then 
the energy transfer is non-local. This is the case when 
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c x n cr© where ct = (m N /M )i?| ~4x 10 -36 cm 2 is 
a critical scattering cross-section. We consider the ADM 
trapped in the Sun as an isothermal gas at temperature 
T x [I], so the luminosity L x carried by the particles is: 

L x (r)= f dr> 4nr' 2 p(r')e x (r') , (6) 
Jo 

where e x (r') cx (T(r') — T x )N x a x ^ is the energy trans- 
ferred to the ADM per second per gram of nuclear matter 
and p{r') is the density in the Sun [Tj. The ADM temper- 
ature T x is fixed by requiring that the energy absorbed 
in the inner region (T(r) > T x ) is equal to that released 
in the outer region (T(r) < T x ), such that L X (R Q ) = 0. 

This approximation overestimates the energy transfer 
by a small factor [55] 122] but is sufficiently accurate for 
the present study. We adopt a simple polytropic model 
for the Sun's temperature T, number density n p and 
gravitational potential V [TJ. The resulting variation of 
the solar luminosity 5L(r) = L x (r) / Lq(t) is shown in 
Figj2] assuming ct x n = 4 x 10~ 39 cm 2 (i.e. 10~ 3 <7 Q ) and 
N x = 2 x lO~ n A from Fig[lJ Note that the luminosity 
scales linearly with both ct x n and N x /Nq. 
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FIG. 2: The radial variation of SL(r) = L x (r)/Lo(r) due to 
ADM of mass 5 GeV, using the approximation of Ref. [T] and 
L Q (r) from the BS05 (OP) Standard Solar Model [24]. 

From the radiative transport equation it follows that a 
small variation of the solar luminosity is equivalent to an 
opposite small variation in the effective radiative opac- 
ity: SL(r) Skj (r) = -K x (r)//t 7 (r) [30] ■ The effect of 

such a localised opacity variation in the region r < O.2i? 
has been studied by a Monte Carlo simulation [31] and 
results in excellent agreement obtained using a linear ap- 
proximation to the solar structure equations [35] . Fig. [2] 
shows that the opacity modification due to a 5 GeV ADM 
with a relative concentration of 10~ n is roughly equiva- 
lent to the effect of a 10% opacity variation. The impact 
of this luminosity variation on neutrino fluxes can be es- 
timated by evaluating SL(r) at the scale height of the 
ADM distribution, 5L(r x ) [TT]. In general, to have an 
observable effect requires o- x nN x /o-qNq > 10~ 14 . 



It is possible through helioseismology to determine e.g. 
the mean variations of the sound speed profile {5c/ c) and 
density (6p/ p) of the Sun, as well as the boundary of the 
convective zone Rcz- In particular Rcz is determined 
to be (0.713 ± 0.001)i? o while the SSM with the revised 
composition [57] predicts values that are too high by up 
to 15(7 |26j . Lowering the opacity in the central region 
of the Sun with ADM also lowers the convective bound- 
ary. The 10% opacity variation shown in Fig. pleads to a 
~ 0.7% reduction in Rcz [35] and thus restores the agree- 
ment with helioseismology. The sound speed and density 
profiles which are presently underestimated in the region 
0.2i? Q < r < Rcz would also be corrected by the opacity 
modification displayed in Fig. [2j 
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FIG. 3: Neutrino producing regions in the solar interior 

The modification of the luminosity profile extends into 
the neutrino producing region as displayed for the SSM 
in Fig. [3] Comparing with Fig. [2] we see that preci- 
sion measurements of different neutrino fluxes may be 
able to test the ADM model and determine its param- 
eters. The ADM mass determines the scale height r x , 
hence the relative modifications of individual neutrino 
fluxes, while the cross-section determines the capture 
rate and thereby the overall modification. Both Monte 
Carlo simulations [31] and the 'linear solar model' [32] 
show that the variation of neutrino fluxes with respect 
to localised opacity changes in the neutrino producing 
region (r < O.2i? ) scales approximately as 6$b ~ 1.55k 
and 5&Be ~ 0.75k. The opacity variation in Fig. [2] 
leads to variations <5$b = — 17%,(5$Be = —6.7% and 
<5$ N = -10%,(5$ o = -14% [32]. Measurements of the 
8 B flux by Super-Kamiokande [33J, SNO [3J] and Borex- 
ino [3S] are precise to 10% while the expectations vary by 
up to 20% depending on whether the old [5S] or the new 
[27] composition is used [26]. For the 7 Be flux, the theo- 
retical uncertainty is 10%, while Borexino aims to make 
a measurement precise to 3% [35]. SNO+ is expected 
to make a first measurement of the pep and CN-cycle 
fluxes [37] . Thus the effects of metallicity and luminosity 
variations can be distinguished in principle. 
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CONCLUSIONS 

Asymmetric dark matter does not annihilate upon cap- 
ture in the Sun and can therefore affect heat transport in 
the solar interior and consequently neutrino fluxes. This 
is particularly true for particles with self-interactions 
which would also explain the paucity of sub-galactic 
structure. We have shown that the presence of such parti- 
cles in the Sun can solve the 'solar composition problem'. 

Intriguingly a 5 GeV 'dark baryon' would naturally 
a) have the required relic abundance if it has an initial 
asymmetry similar to that of baryons, b) have a self- 
interaction cross-section of the right order to explain sub- 
galactic structure, c) modify the deep interior of the Sun, 
restoring agreement between the standard solar model 
and helioseismology, and d) be consistent with recent 
hints of signals in direct detection experiments [T7J [05] . 
Such a 5 GeV ADM particle would lower the solar neu- 
trino fluxes which ought to be testable by the Borexino 
and (forthcoming) SNO+ experiments. 

Note added: After this paper was submitted to arXiv 
(1003.4505), another study appeared [35] with similar 
findings concerning the effect of ADM on helioseismol- 
ogy and neutrino fluxes. However a second such study 
[40] finds negligible effects using a solar model simulation. 
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